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Abstract 

We calculate correlation functions in matrix models modified by trace-squared 
terms. First we study scaling operators in modified one-matrix models and find that 
their correlation functions satisfy modified Virasoro constraints. Then we turn to 
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tions provide additional support for the idea that the modified matrix models contain 
operators with the negative branch of gravitational dressing. 
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1 Introduction 



In recent literature some effort has been devoted to studying the large N matrix models 
modified by trace-squared terms []. In the one- matrix models, for instance, one may add 
terms like g(Tr $ 4 ) 2 to the conventional matrix potential of the form NTt V(<&). 

Discretized random surfaces appear in the Feynman graph expansion of the matrix mod- 
els, and the vertex corresponding to the trace-squared term may be thought of as an iden- 
tification between a pair of ordinary plaquettes. If the two plaquettes belong to otherwise 
disconnected random surfaces, then this identification introduces a tiny "neck" into the ge- 
ometry; if they belong to the same connected component, then the identification attaches a 
microscopic handle. These two types of effects are familiar consequences of the Euclidean 
wormholes, which were widely explored in the four-dimensional quantum gravity. Thus, 
trace-squared terms are a convenient trick for adding wormholes to the random surface 
geometry, with the correct combinatorics emerging automatically. Actually, wormholes are 
known to be abundant already in conventional random surface theories with no trace-squared 
terms []. Therefore, turning on the coupling g does not introduce any new types of geometry 
into the path integral, but simply increases the weight of microscopic wormhole configura- 
tions. For example, in the sum over surfaces of genus zero there are geometries corresponding 
to trees of smooth spheres glued pairwise at single plaquettes, and increasing g enhances their 
weight. It is not surprising, therefore, that a small increase in g does not change the universal 
properties of the model. If, however, g is fine-tuned to a finite positive value g t , then the 
universality class of the large area behavior changes. For pure gravity (c = 0) the string 
susceptibility exponent jumps from —1/2 for g < g t to 1/3 for g — g t . This is the simplest 
example of a matrix model where new critical behavior occurs due to fine-tuned wormhole 
weights. 

Further work has revealed that, more generally, as the trace-squared coupling is increased 
to a critical value g t , the string susceptibility exponent jumps from some negative value 7, 
found in a conventional matrix model, to a positive value 

< L1 > 

Essentially equivalent results have been obtained without using matrix models, on the basis 
of direct combinatorial analysis [] . For a long time the positive values of string susceptibility 
exponent seemed very puzzling. Recently, however, a simple continuum explanation of these 
critical behaviors was proposed in ref. []. 

For all the conventional matrix models describing (p, q) minimal models coupled to grav- 
ity, the correct scaling follows from the Liouville action of the form 



i J d 2 <7 ((d„4>f - QR4> + (O min e°+*) 



° + = M^- C+ 24,v "' - ^) = J1± Mr • (L2) 

where O m i n is the matter primary field of the lowest dimension. A simple calculation reveals 
that the string susceptibility exponent is given by 

7 = 2 + ^, (1.3) 

a + 



1 



where Q = y In ref. [] it was argued that the effect of fine-tuning the touching in- 

teraction is to replace the positively dressed Liouville potential by the negatively dressed 
one, 



p + q + 1 



1 - c + 24fr min + \/25^ =- - L- ■ (1.4) 



2^3 V ' mm ' V J 

Now the string susceptibility exponent is found to be 

7 = 2 + ^ = ^-, (1.5) 
a- 7 — 1 

in agreement with the matrix model results. The intuitive reason for the change of the branch 
of gravitational dressing is that the microscopic wormholes alter the ultraviolet (large 0) 
behavior of the theory. A priori, there are two independent solutions to the string equations of 
motion, corresponding to the two choices of dressing. The linear combination that appears in 
the theory is selected by the boundary conditions, and by fine-tuning the boundary conditions 
at large <fi we have changed the string background. 

A more complete understanding of the modified matrix models was recently found in 
ref. []. As the trace-squared coupling is set to g t one finds the following non-perturbative 
relation _ 

- _ roc _ 

e F(t;tl,t2,...,t n ) = / fa e tt+F{t-,tx,t2,...,tn) ( h 5) 



where F and F are the universal parts of the modified and conventional free energies, re- 
spectively. As indicated in eqs. ( |1.2|) and (|1.4j) , t is the lowest dimension coupling in the 
conventional action and i is the corresponding coupling in the modified action; ti are coupling 
constants corresponding to other operators. Evaluating the integral in (|1.6| ) in saddle-point 
expansion, we may generate the genus expansion of F in terms of the known genus expan- 
sion of F. This gives a concise prescription for calculating universal correlation functions in 
modified matrix models, which will be used extensively in this paper. 

A further result of ref. [] is that, by changing the type of trace-squared terms, it is 
possible to introduce integrations over other couplings U. For example, if t\ is the coupling 
constant for a matrix model scaling operator 0\, then by adding g\0\ to the matrix action 
and fine-tuning gi, we obtain a model where 

roc 

e F(t;ti,t2,-,t„) _ / C ^ ie titi+-F(t;tl,t2,...,tn) _ 
J — oo 

The dependence of F on t\ indicates that the gravitational dimension of 0\ has changed 
from its conventional value d to 

J=7-d (1.8) 

(the string susceptibility exponent remains unchanged). Remarkably, this change of dimen- 
sion is again reproduced in Liouville theory by a mere change of the branch of gravitational 
dressing. Namely, if the operator is dressed by er^ then 

d — 1 — -^t, d = 1 — — . 
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The relation ( |1.8|) follows from j3 + + /3_ = — Q. 

By applying integral transformations ( |1.6[) and ( |1.7|) sequentially to any subset of cou- 
pling constants, we can construct a matrix model where the scaling dimensions of all the 
corresponding operators correspond to picking the negative branch of Liouville dressing. It 
appears that we have found a whole new class of matrix models which serve as exact solu- 
tions of a new class of Liouville theories, those involving some number of negatively dressed 
operators. The purpose of this paper is to use the new matrix models to extract some in- 
formation about the negatively dressed operators and to compare it, if possible, with direct 
continuum arguments. 

In section 2 we analyze the correlation functions of macroscopic loop operators w(£) in 
modified one-matrix models. Using eq. (|1.6|) to implement the change of dressing of O m j n , 
we calculate the loop correlators on a sphere and find agreement with the results of refs. []. 
Correlation functions of the scaling operators a n may be read off from the expansion of w(£) 
in powers of t. In addition to the basic operator o§ which by definition has gravitational 
dimension 0, we find operators a n , n > 0, of dimensions j(n + 1). Remarkably, there is no 
operator of dimension 7 which in the Liouville language would correspond to O m i n e a+< ^. This 
is not a coincidence: operator O m i n e a -^ already appears in the theory, and a simultaneous 
appearance of O min e a+<?i would signify a doubling of the spectrum that is unacceptable on 
general grounds. [] 

The dependence of F on t n , the coupling constants for a n , can be found from the gen- 
eral formula ( |1.6| ) (an explicit derivation of this fact for this particular system appears in 
Appendix A). From the fact that e F obeys the Virasoro constraints, it follows that e F obeys 
modified Virasoro constraints, obtained from the conventional ones by replacing 

We discuss the Virasoro constraints and the recursion relations they imply among the cor- 
relation functions in section 3. 

In section 4 we check formula ( |1.7| ) for transformation of scaling operators in the context 
of one-matrix models. For this formula to apply to gravitational descendants, a class of 
leading analytic terms must vanish in their conventional correlation functions. We check 
this vanishing for some specific examples. 

In section 5 we use eq. ( |1.6| ) and formulae for the spherical two- and three- point cor- 
relators of positively dressed c < 1 order parameters to calculate correlators involving the 
negatively dressed order parameters. We find simple results which agree with a straight- 
forward analytic continuation in the Goulian-Li formulae []. As a further check, we carry 
out similar analysis for the c = 1 theory up to the four-point function which is, fortunately, 
known explicitly. 



*As we remarked before, in any theory we expect the boundary conditions to single out unique gravita- 
tional dressing. 
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2 Scaling Operators and Macroscopic Loops 



In this section we investigate the correlation functions of scaling operators in the modified 
multicritical one-matrix models [1, 



X)$ e - Ar ( Tryfc (*)+( C2 - A ) Tr * -2Tv( Tr * ) ). (2.1) 

The critical potential of the kth model with g = is 

k 

y fe ($) = 5>$ 2 * , (2.2) 

where c« have been determined in ref. []. 

The correlation functions of the scaling operators for g = have been studied in ref. []. 
The scaling operators are given by linear combinations of traces of powers of $, 

n+1 

cr n = ]T^ n) Tr$ 2 \ n = 0,1,2,... . (2.3) 

i=i 

The coefficients <?} are chosen so that the correlation functions of cr n 's on the sphere scale 
as []: 

{a x ... a p ) ~ i^T+EnC**- 1 ) , (2.4) 

where t ~ (c2 — A)iV 2// ( 2-7 ), and d n = — n7 is the gravitational dimension of the operator cr n . 

This analysis can be extended to the g ^ case. For the fine-tuned g = g t , in addition 
to the dimension operator, &o = Tr$ 2 , one finds a set of "modified" scaling operators 

n+1 

a n = J29i n)r R® 2i , n=l,2,3,..., (2.5) 

i=i 

with gravitational dimensions c? n = (n+1) 7. The correlation functions of these new operators 
on the sphere are given by 

ff p ) ~t 2 ^+E„(^n-l) ) ( 26) 

where t ~ (A c - A)A^ 2 /( 2 -^). There are two interesting features of these results we want to 
emphasize. First, the coefficients g^ in (|2.5| ) are equal to the coefficients gf 1 ^ 1 in ( |2.3| ) for 
all i except i = 2. A general expression for g^ in the fcth multicritical model is given in 
Appendix A. Second, there is no operator of dimension 7. This provides additional argument 
in favor of Liouville interpretation of the modified matrix models []. In Liouville theory such 
an operator would correspond to O m i n e a+< ^ which is not acceptable. Namely, by choosing 
the negative branch of dressing we already have operator O m i n e a_ ^ in the theory, and the 
existence of the operator O ra \ n e a+ ^ would mean a doubling of the spectrum. 

In Appendix A we carefully derive eq. ( |1.6[ ) which generalizes the relation between the 
universal parts of modified and conventional free energies in the presence of perturbations of 
the multicritical potentials. This result allows us to calculate arbitrary correlation function 
of a n 's for any genus once the corresponding correlation functions of cr n 's are known. 
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It is actually more convenient to consider the macroscopic loop operator, w(£), which 
contains the complete information about correlation functions of scaling operators cr n 's in a 
compact form, 

°o pn+1/2 

Let us assume that tj's in eq. ( |l.(j| ) are couplings to macroscopic loops w(£i), and denote 
them by r^. Then we have 



and 



(w(£x)...w(£ p )) 
{w{£ x )...w{£ p )) 



dri . . . dr p 



F(t- n 



F{i;n) . 



(2.8) 



(2.9) 



dr\ . . . dr v 

We now show how correlation functions in the modified models (|2.9| ) can be systematically 
calculated in terms of those in the conventional c < 1 models ( |2.8j ). 

Let us first consider the correlators on the sphere. In this case eq. (|1.6|) can be written 
as Legendre transform 

F(t;r i )=tt + F(t;r i ) , (2.10) 
where the r.h.s. is evaluated at the saddle-point t s (t ), given by the solution of equation 

dF(t; ri ) 



dt 



(2.H) 



By simply taking derivatives of fl2.10|) , and using the chain rule, one finds for the one-loop, 



(w(£))o(i) = (w(£))o 



t=t a (t) 



two-loops, 



(w(£ 1 )w{£ 2 )) {i) = {{w(£ 1 )w(£ 2 ))o- 



t=t(t ) 



three-loops, 



(w{£i)w{£ 2 )w{£z))o{i) = (w(e 1 )w(£ 2 )w(e 3 )) 



(w(e 1 )P)o(Pw(£ 2 ))o 
(PP)o 

(w{£ x )w{£ 2 )P) (Pw{£z)) 



(2.12) 



(2.13) 



(PP)o 



(w(£ 1 )w(£ 3 )P) (Pw(£ 2 ))o (w(£ 2 )w(£ 3 )P)o(Pw(£ 1 )) 



(PP)o 



(PP)i 



_ (w(£ 1 )PP)q(Pw(£ 2 MPw(£ 3 )) (w(e 2 )P P)o(Pw(e 3 )) (Pw(£ 1 ))o 

/ I 7-1 7 "i\ \ O 



+ 



((^)o) 2 
(w(£ 3 )PP) (Pw(£ 1 )) (Pw(£ 2 )) 
((PP)o) 2 

(w(£ 1 )P)o(w(£ 2 )P}o(w(£ 3 )P) {PPP} 
((PP)o) 3 



((PP)o) 



(2.14) 



t=t(t ) 
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and so on. In the above expressions subscript denotes the genus of the surface and P is 
the puncture operator. Using the known results for the loop correlators when g = [] one 
easily checks that our results ( |2.12| ), ( |2.13| ) and ( [2.14| ) agree with the direct calculations of 
loop correlators in refs. [] . 

For example, for the k = 2 model on has explicitly: 



««)> = i(i+«- 1/3 k a """ 



£5/2 



2 



2 



(w(4M4)™(4)> = 7 -7=^ (4W 1/3 + (44 + 44 + 44)r 2/3 

4 V444 

+ (4 + 4 + 4)* -1 + * " 4/3 ) e-^ +/a+ *») f 1/3 . (2.15) 

Expanding (w(£)) in powers of i, one finds that the power y/i is missing, which is related to 
the absence of the dimension 7 operator. Similarly, in the expansion of (w(4)^(4)) there 
are no terms with \fl\ or \/4, etc. Thus, instead of (|2.7|) we may write 

1 OO £71+1/2 

wit) = —a + Y — — a n . (2.16) 

This reflects the fact that the new puncture operator comes from analytic terms in the 
conventional model, while the conventional puncture operator becomes analytic. 

Note that the relation Q2.10 ) between loop correlators in modified and conventional matrix 
models is the same as the relation between the generating functions of connected and one- 
particle irreducible Green's functions. This analogy suggests a simple diagrammatic relation 
between the two sets of correlators. 

Let us now extend the above analysis to surfaces of arbitrary genus. The saddle-point ex- 
pansion of eq. ( |1.6p generates the complete genus expansion. It is, however, more convenient 
to formulate a set of graphical rules ("Feynman rules") which allow us to relate correla- 
tion functions in two models in a simple and geometrically transparent way. Expanding the 
exponent in the integrand of eq. (|1.6|) around the saddle point, t = t s , we have, 



e p = [ dtexp \ t A + l A t + F(t s ) + -(P)At+ V^-(P n )(Atf 1, (2.17) 

J I K K K n! K I 



where 



A* = *-*,, -(P n ) =F^(t s ) , 

K 

and we have exhibited explicit dependence on the string coupling constant, k. The saddle 
point is determined by 

t+(P) = 0, (2.18) 

where (P) is the exact one-point correlation function of the puncture operator. Since we are 
eventually interested in the expansion of e F in powers of k, it is most convenient to solve 
eq. ( |2.18| ) on the sphere, 

t+(P)o = 0, (2.19) 
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<pp> f 




K g_1 <P> a . g=l,2,3, 





V^, ... ) — K g - ! <pp>„, g =i,2,3, 




K g-1 / p n\ , §=0,1,2,3,- 
K {F } s n=3,4,5,... 



Figure 1: Graphical rules for constructing correlation functions in the modified matrix models 
in terms of those in the conventional models. 



and introduce explicit tadpole terms for surfaces of genus g > 1. Similarly, the propagator 
which we read off the quadratic part in eq. ( [2.171 ), — K / (PP), depends on k in a complicated 

way, 

_ ^_ = 1 = _ V :- ^ PP)g 

(PP) E g > ^(PP) 9 (PP)o\ g (PP)o 

(PP) gi (PP) 

Therefore, we choose as a propagator — k/ (PP)q, and take higher-genus mass insertion terms 
as vertices. The complete list of graphical rules is shown in fig. |l[ 

As an example, we show the one-loop correlator on the torus (g = 1) in fig. 0. The 
corresponding analytical expression reads: 

(^.-^-fi^-i^.ite^Qs. „ 21 > 
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w(l)) = 

' i 




+ 




+ 




+ 




Figure 2: One-loop correlator on the torus. 

The reduction formulae do not work for correlation functions involving puncture operator. 
Instead, one uses the chain rule which follows from eq. Q2.19 ): 



v ' dt n N ' \dtdt s 




with particular cases 



and 



d d 

(P) = F = — (it s + F (t s )) =t s = t 1 ' 



(PP)i 



1 



(PP)o 

A higher genus example is the one-punctured torus: 

(P) =^F - {P)l I 1 {PPP)o 

[ n at 1 (pp) 2 (pp) 



(2.22) 

(2.23) 
(2.24) 

(2.25) 



3 Virasoro Constraints and Recursion Relations 

One of the most important mathematical properties of the c < 1 matrix models is the 
existence of Virasoro constraints on the partition function (or appropriate generalizations for 
multi-matrix models). These constraints are equivalent to the loop equations and are related 
to the integrability property of these models, both on the lattice and in the continuum. When 
expressed in terms of correlators of local scaling operators, they take the form of recursion 
relations identical to those defining topological two-dimensional gravity. A nice geometrical 
picture arises at the multicritical points with 7& = —1/k: all correlators of scaling operators 
cr m with m > k — 2 reduce to contact terms at the boundaries of moduli space and can 
be solved in terms of the correlators of o m with m < k — 1. This defines the notion of 
gravitational primaries and descendants, a gravitational version of relevant and irrelevant 
perturbations of the critical point. 

The Virasoro constraints in the continuum take the form [] (defining t = t): 



L n Z{t ,{t}) = , n>-\ 



(3.1) 
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where {t} = ti,t 2 , ■ ■ ., and Z = e F is the disconnected sum over continuous surfaces.^] The 
differential operators L n are given by 

^o= e r + 2j tm ar + T6 ' 

m>0 V Z/ Lam ±U 

L « = E fm+-)t m ^^ + -E ^ , (3-2) 

where we made explicit the dependence on k, the string loop expansion parameter. These 
operators satisfy a centerless Virasoro algebra 

[L n , L m ] = (n - m)L n+m . (3.3) 

In the modified matrix models, the partition function (and in general any disconnected 
correlator) is defined in terms of the corresponding object in the standard matrix model by 
the Laplace transform 

Z(t , {*}) = J dt e toI °/ K Z(t , {£}) = C[Z] . (3.4) 

This formula makes sense as the saddle point or genus expansion and, as was discussed 
in the previous section, any correlation function of loops or scaling operators in models of 
"touching" surfaces may be decomposed into sums of products of correlators of conventional 
models. It is then clear that the conventional recursion relations imply certain modified 
recursion relations in the new models, perhaps with a similar topological interpretation. In 
fact, such relations follow immediately from ( |3.4[ ) if we define a set of modified Virasoro 
operators by the operator identity 

L n C = CL n , (3-5) 

thus satisfying 

L n Z{i ,{t}) = , n>-l. (3.6) 
The new operators are related to those in ( p.2|) by the transformations 

*o -» K 7jr > ~W -> • ( 3 - 7 ) 

One can explicitly check that the L n operators so defined satisfy the same Virasoro algebra 
as L n . A simple proof of this fact follows from the abstract Fock representation introduced 
in ref. [] . The Virasoro operators (|3.2|) can be interpreted in terms of the energy- momentum 
tensor of a twisted scalar field on the circle, with mode expansion 

^)=E«n + I^ n - 3/2 , (3-8) 



2 The continuum limit of matrix models with even potentials induces a doubling of the degrees of freedom, 
so that the Virasoro constraints act on the square root of the matrix model partition function in that case. 
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under the following substitutions for n > 0: 

a-n-f = ^ (n + i) t„ , a n+| = • (3.9) 

It is easy to see that the Fock space representantion of the modified Virasoro operators is 
linearly related to the previous one, a n +i/2 = «n+i/2 for n / 0,1 and 0:1/2 = — 2o_i/2, 
a_i/2 = |«i/2- In fact, it is a Bogoliubov transformation, since it conserves the canonical 
algebra 

[a n+ i,a m _i] = [a n+ i,a m _i] = (n + -) <5 n+ m,o (3.10) 

from which one derives the Virasoro algebra ( |3.3| ) for the L n operators. 

Recursion relations are easily obtained starting from the general identity: 

n|-(e-^»e p ) = 0, (3.11) 

expanded in powers of the string coupling at the kth multicritical point, t n = —5 ni k/(2k + l), 
n > 1, with string susceptibility % = l/(k+ 1). Let us consider for simplicity the case in 
which there is no puncture operator a = P in the set S, and use the notation &g = Yl s aj. 
Neglecting some analytic terms in the couplings, the L_i or puncture equation takes the 
form: 

((Tk-i*s) a = £(2j + + \(PP<ts) g -i + ^ £ £ ^' v (/Vv), /: , • (3.12) 

jeS 4 8 XUY 91 +92=3 

The first term on the r.h.s. represents operator contact terms and is identical to the conven- 
tional counterpart. However, in the modified model we have additional factorization terms 
from the contribution of the puncture operator at the boundaries of moduli space, where a 
genus g surface degenerates into two g\ + g 2 = g surfaces or a g — 1 surface by pinching a 
handle. 

The Lq or dilaton equation is identical to the conventional one, up to a sign flip of the 
bulk term, proportional to to'- 

(ws) B = E(2j + m^S-i) g - i (Pcrs) 9 • (3.13) 

It is interesting that this equation is insensitive to the puncture factorization. The higher 
n > 1 equations take the following form: 

(a k +n&s) g = £(2j + l)(a j+n a S -j)g - 2t (l - 5 n) i)(a n -ias) g 

n-l 

+ (P&n&s)g-1 + £ (&m-l&n-m&s)g-l 
m=2 

• (3-14) 

Z Xuy3i+32=3\ m=2 / 
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Again, the operator contact term as well as the factorization terms not involving the puncture 
operator are identical to the conventional case. The factorization for the puncture is different 
though: it comes with an extra factor of 1/2, and the conjugated operator at the degeneration 
neck is a n instead of a n _i. The bulk term is also slightly different. 

These equations are still recursion relations for the gravitational descendants, in the usual 
sense, and can be derived from the corresponding relations in the conventional models using 
the diagrammatic method explained in the previous section. This represents a non-trivial 
check of both the recursion relations and the diagrammatic rules for touching surfaces, due 
to the complicated pattern of cancellations involved. For this reason it is perhaps surprising 
that the final result is so similar to the conventional case, with only some modifications 
in the way the new puncture operator enters in factorization diagrams and bulk terms. We 
interpret this fact as evidence for an underlying topological field theory description, involving 
only slight changes from the usual one, perhaps similar to the change of dressing branches, 
advocated for the Liouville model. Presumably, the role of the pure topological point, k — 1, 
is played here by the pure polymer phase with 71 = 1/2. 

To conclude this section, we remark that one can also write integrated forms of the L_i 
and L equations acting on loops, measuring boundary lengths and overall dilatations. To 
this end, we simply insert the loop operator 

1 d £«+i/2 q 

^ e) = 7tw a + ^f^rm^ (3 - 15) 



on the left hand side of eq. (|3.11|) 



We find some differences with respect to the conventional models. The dilaton operator, 
which at the kth critical point is a^, satisfies 

ifo*M> = (-i«,£+4)<«(0>. (3-16) 

Comparing with the conventional models, the sign of the bulk term is flipped. In the modified 
models we find additional non-local terms associated with the new puncture contributions 
to the L_i constraint, as well as an inhomogeneus term, proportional to the one-puncture 
function, 

^(a k ^w{£)) = £(w(£)) - y/i(P) + i (k(PPw{1)) + 2(P)(Pw(£))) , (3.17) 

where we have neglected some analytic terms. Thus, it seems that a local boundary operator 
is absent in the modified models. (<Tk-i acts as a boundary operator in the conventional 
models due to the relation ^((Jk-iw(£)) = £{w (£)).) There are interesting generalizations of 
the above exercise. For example, one may consider H^-constraints for modified multi-matrix 
models. Also, other operators can be tuned as wormhole sources. In the most general case 
one has a multiple Laplace transform and still a Bogoliubov transformation in the Fock space 
representation. 



11 



4 Changing Gravitational Dimensions 



In the previous sections we studied modified one-matrix models with the simplest type of 
trace squared term, (Tr$ 4 ) 2 . We showed that the effect of this term is to change the branch 
of gravitational dressing of the lowest dimension operator in the theory, which corresponds 
to coupling constant t. In this section we check that more complicated trace-squared terms 
alter the gravitational dimensions of other scaling operators. In particular, we check in detail 
that formula ( |1.7| ), first derived in ref. [], applies to gravitational descendants. 

Although our approach is general, we focus for simplicity on the k = 2 one-matrix model, 
whose matrix potential is 

S ($) = NTt (^$ 2 - A$ 4 ) . (4.1) 

The continuum limit is achieved as A = A c — A — > 0. It is convenient to introduce the first 
scaling operator of the form 



a 1 = NTt (V - ^$ 4 ) - N 2 (l - 16A - 2304A 2 ) . (4.2) 

This operator has gravitational dimension d — 1/2, and its connected genus zero correlation 
functions are given by 

(ax) = iV 2 0(A 5/2 ) , 

(<7 lt7l ) = N 2 (^ - 1024V/3A 3 / 2 + 0(A 2 )) , 
(o-xaKJi) = iV 2 (1536A + 0(A 3/2 )) , 

(<7i<7i<7i<7i) = N 2 (-7Q8V3A 1/2 + 0(A)) , etc. (4.3) 

The purpose of the ^-independent term in o\ is to remove a non-universal analytic term 
from the one-point function. Apart from this it has no effect. 
Now we consider a modified model with the action 

So($)-^i-^(<ti) 2 , (4.4) 

where we have introduced coupling constant T\ in order to study correlation functions of o\. 
Its partition function may be written as 

Z ~ l°° dve'^f [ V$e- So+{v+Tl)ai . (4.5) 

Defining a shifted variable u = v + T\, we perform the matrix integral first and reduce the 
modified free energy to 



F = F (t) +log J due fiu) 



N 2 4 
f(u) = - — (u 2 - 2un + t 2 ) + -N 2 u 2 + F^h) , (4.6) 
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where the scaling variables t and t± are defined through 

t ~ AiV 2 /( 2 -^ , t X = u iV2(l-d)/(2- 7 ) (4 ?) 

In this specific case 7 = — 1/2 and d = 1/2. Of crucial importance is the fact the F\ depends 
only on the scaling variables and has the form 

00 

F 1 (t,t 1 ) = E^ (5 " n)/2 > ( 4 -8) 

n=2 

where a n are constants (the first three of them may be read off eq. fl4.3j )). This follows from 
miraculous vanishing of certain potentially harmful non- universal terms in eq. (|4.3|) . For 
example, had the three-point function of o\ started with a non-universal term of order A , 
eq. ( |4.8| ) would not be valid! 

If we now set g = 3/8 and introduce the scaling variable 

t, = 8 TliV 2(l-<i)/(2- 7 ) ; (4 g) 

3 

we arrive at the following expression for the universal part of the modified free energy, 

/oo _ 
dti e titi+F(MO j (4.10) 
-00 

Here F(t,ti) is the universal part of the conventional sum over surfaces. Thus, we find that 
the gravitational dimension of o\ has changed from d = 1/2 to d — j — d — —1. This 
provides a counterexample to the claim of ref. [] that eq. ( |1.7| ) applies only to operators with 
d<(l + 7)/3. 

While in ref. [] eq. (1/7 ) was derived for gravitational primary fields, we have just checked 
it for a gravitational descendant. Our example shows that these operators have to be defined 
in such a way that their correlation functions do not contain certain non-universal contribu- 
tions. We believe that this is in fact always possible. As further evidence, we present results 
for the dilaton operator, 

a 2 = NTi ( 2$ 2 - V + -$ 6 ") - N 2 (- - —A] . (4.11) 
V 3 60 / \5 5 / 

This operator has gravitational dimension d = 1, and its connected genus zero correlation 
functions are given by 

(a 2 ) ~ iV 2 A 5 / 2 , 

(^ 2 ) = iV 2 (| + 0(A 5 / 2 )), 

(cr 2 a 2 a 2 ) ~ iV 2 A 5 / 2 , 

\o 2 o 2 o 2 o 2 ) ~ iV 2 A 5/2 , etc. (4.12) 

Once again, the unwanted non-universal terms vanish! Repeating the steps carried out for 
01, we may now establish the validity of ( |1 . 7| ) for a 2 . This relation, and the change of 
gravitational dimension it implies, appear to be completely general. 
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5 Correlators of Dressed Primary Fields 



In the previous section we studied correlation functions of gravitational descendants in mod- 
ified one-matrix models. In this section we turn to gravitationally dressed primary fields. 
Although our approach is general, we mostly discuss some simple special cases: genus zero 
two- and three-point functions of the order parameter fields in unitary minimal models cou- 
pled to gravity, as well as the four-point function for c = 1. We find that the correlators 
in modified matrix models agree with those of negatively dressed operators in Liouville the- 
ory, provided that the latter are obtained with the simplest plausible analytic continuation 
prescription. Although a far better understanding of the Liouville theory calculations is 
desirable, we feel that our findings suggest a general pattern for interpreting the modified 
matrix models in terms of the negatively dressed operators. 

The starting point for our calculations is the relation ( |1.7p between the genus zero free 
energy F of the modified matrix model and that of the original matrix model, F: 

dF 

F(t 1 ,i2,...) = ^2t i i i + F(t 1 ,h,...), — = -tj , (5.1) 

where ti and U represent the coupling constants by which the model is perturbed (it is 
implicit that F and F also depend on the basic coupling constant t, which appears in the 
Liouville action). Since ti correspond to dressed primary fields, F has the expansion 

F(ti) = l -^tl + l -c lJk nn k + ±d m tHH k t l + ... (5.2) 



where Aj, Cjj&, diju are the ordinary two-, three- and four-point functions. F in (|5.1| ) can be 



evaluated order by order in the U by expanding UU + F(U) around the saddle point. 

Let us begin by discussing the two- and three-point functions. First assume that there 
is only one parameter £j. Then the Legendre transform fixes £j at 

ti0 = ~ A~ U ~ 2Af ** + --" (5 - 3) 
The genus zero part of F is just the value of tjij + F(ti) at t i0 . One finds in this case 

F(ti) = -A,t7 + \c m if + ... (5.4) 
2 o 

with modified two- and three-point functions 



\ — ^ j c m — ^3 . (5.5) 



Next, assume that there are two parameters ti, tf. and one wants to modify one of them, i.e., 
one wants to find F(ti,tj.) to cubic order. Now the saddle point is at 



HQ 




C %t k ) ti + . . . (5.6) 
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In this case, 

t k ) = \^tl + \-A k t\ + \c Tlk tHH k + \c- lkk f t k t k + . . . (5.7) 
2 2 o o 

with 

i = T~~ j Cuk = ~r a 2 ? Cjfcfc = — . l^-"J 

ZAj ZAj ZAj 

The generalization is straightforward: each time we Legendre transform from ti to U the 
two- and three-point functions change according to 

A« > Aj = — — , Cij k ► Cij k = . (5-9) 

The absolute values of the correlators are normalization dependent, so it is useful to define 
the normalization independent quantity as in ref. []: 

c 2 

X ijk = % J\ F , (5.10) 
l3k AiAjA k v ; 



where F is the free energy evaluated at ti = 0. One easily sees from eq. (|5.9| ) that X has 
a simple property: it just switches sign each time one of the external operators is modified, 
i.e. 

Xijk = X{j k , Xij k = X{j k , Xy k = Xij k . (5.11) 

Let us compare this behavior with what one would expect from Liouville theory, if the 
modified matrix model operators (5j were identified with the negatively dressed operators. 
Consider the gth minimal model with central charge 

6 



q(q + l) ' 

coupled to gravity. The gravitationally dressed operators on the diagonal of the Kac table 
are 

Oi = A it n ^ U , O t = Aun , (5-12) 

with 

,± Q 



2 ' ' J2q(q + 1) 



We note that switching from the (3f to the f3~ dressing corresponds to continuing — > —r^ 
The dressed unity e a+ ^ corresponds to the puncture operator with = 1. As shown in 
section 2, if we Legendre transform with respect to the cosmological constant t, as in eq. 
( p,10| ), we find a new theory with Liouville potential e a ~^. 

First, note that the scaling behavior of the modified correlators (|5.5|) and ( |5.8|) is that of 
operators with the /3~ dressing. This has already been noted in ref. [] and in section 4. Next, 
let us compare the coefficients. Consider the quantities Xij k , defined in ( |5.10|) . They can be 
derived from the fact that the three-point function of the properly normalized operators can 
be written as [] 

Q < 3 t+r ) t+' 3 i 



{Oi0 3 O k )=t <*+ . (5.13) 
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Since insertions of the puncture operator P are produced by —d t , we can infer from ( p. 13 ) 
and from the scaling behavior the two-point function and the partition function: 



-2^—1 , „ _ , OL+ _-2-_ 2 3L 



and 



d t (O k O k ) = (PO k O k ) = t *+ (O fe O fe ) = n + + r° + " "+ , (5.14) 



a 2 F = (PP) = — ^ — t a + F = - 7TTT) -t a + . (5.15) 

1 ^ 1 Q + 2a + Q- -2. + 1 (-2-4-2 v ; 

We thus obtain 

(OiOjOtfF _ (Q + 2f3t)(Q + 2#)(Q + 2/3+) 

A ^' fc (0A><cW<<m> g(g + « + )(g + 2« + ) ^.ioj 

This is the result of Goulian and Li []. In order to calculate Xy k we assume that it is given by 
the above equation with (3 + replaced by /3~. This procedure is consistent with all the tricks 
used in the Liouville theory calculations []. The result is remarkably simple: X^ k = —Xij k , 
in agreement with ( |5 . 1 1| ) . In general, Xij k simply switches sign each time one of the Liouville 
dressings is modified from j3 + to (3~ . It appears that we have found an explanation for the 
changes of sign caused by Legendre transform ( |5.1| ). Remarkably, they have the same effect 
on correlation functions as changes in the sign of the Liouville energies. 

Equally remarkable is the effect on (|5.16| ) of Legendre transforming with respect to the 
cosmological constant t only, as in eq. (|2.10|) . In this modified matrix model we find 

{Cty£ , (5.17) 



Q_(Q_ + 2 ) 

where C = (— + 1)(— + 2). Correlation functions of order parameters other than the 
puncture become "one-puncture irreducible" and may be computed using the same rules as 
in eqs. ( |2. 12| ) — ( p. 14[ ) . The modifications are trivial because one-point functions vanish, and 
we arrive at 

{O k O k ) = - j + (Ct) «- a «- , 



Q + 28 



k 



{0,0,0,) = (Ci) ~- ~- . (5.18) 

This leads to 

(OiOjOtfF (Q + 2(3f)(Q + 2(3f){Q + 2/3+) 



Y 

(Q + 2/3+)(Q + 2/3+)(Q + 2/^ 



(0^(0,0,) (O k O k ) Qa + (Q + 2a + ) 

+ y.ft+VO 4- 9.13+) 

(5.19) 



Q(Q + a_)(Q + 2a_) 

Thus, the sole effect of the Legendre transform with respect to t is to replace a + by a_ in 
eq. ( [5.161 )! This precisely agrees with the idea that the dressing of the Liouville potential 
has changed from e a+ ^ to e a ~^. 
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To generalize our observations to the four-point function, it is useful to note that the 
relation between F(ti) and F(ii) is the same as that between W(j), the generating function 
of the connected diagrams, and S(<f>), the effective action, 



dj e 



j<t>+W(j) 



(5.20) 



At tree level S(4>) is identical to the generating function T(0) of one-particle irreducible 
diagrams. Based on this analogy, we can graphically represent the relation between modi- 
fied correlators (black circles) and original correlators (white circles) in terms of Feynman 
diagrams. Let us illustrate this with some examples: 



(A) 



J 




. \ _ 

■k = ry^—k 

j 



(B) 



z I 
j k 



X 



z X 
m y 




(C) 




(D) 



O + O— •— O + crossing terms 

j k j k 



(E) 



One indeed recognizes formulae ( |5.5| ) in the first two lines (A), (B). Example (C) refers to 
the four-point function Legendre transformed only with respect to a coupling corresponding 
to an operator appearing in the intermediate state (the couplings corresponding to external 
legs remain untouched). This changes the value of the four-point function according to 

dijkl * dijki Cjj m ~ c m ^;. (5.21) 
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The implications of this are quite significant. In the language of Liouville theory, we have 
not modified the action, nor have we modified any of the operators entering the four-point 
function. Nevertheless, the value of the correlator changed. A similar effect was observed in 
[] in the course of calculating the genus one free enegy: it changed in response to changing 
the dressing of any operator, even though the Liouville action remained untouched. This 
suggests a remarkable subtlety in the continuum Liouville calculations, related perhaps to 
contributions from boundaries of moduli spaces. 

In figure (D) we show a simpler example where a four-point function is Legendre trans- 
formed only with respect to a coupling corresponding to an external state, so that 

dijki = —-r-dijki ■ (5.22) 

Aj 

Finally, in example (E) we demonstrate the most general transformation of the four-point 
function under the Legendre transform of some set of coupling constants, 



d ^ = AiAjAftA ~ {? ^^ m h CmH ~kk + Cr ° SSing termS ^ • (5 - 23) 



The sum over m runs over all intermediate operators whose coupling constants are Legendre 
transformed. 

The above discussion suggests that it is useful to introduce the normalization independent 
quantity 

X *W = a a ~a a I dijki ~ \ \ ^2c ijm -^—c mkl + crossing terms ) 1 F . (5.24) 
AiAjAjfcAj [ 2 \ m A m J J 



(Note the factor 1/2 in comparing with Q5.23 ) and (|5.21 ).) Here, the sum over m runs 



over all operators. Using (5.9p, ( |5.21| ) and ( |5.22| ) we establish that X is invariant under the 
Legendre transform with respect to any intermediate operator, but switches sign each time 
an external operator is transformed. 

Let us compare this with the behavior of correlation functions in c = 1 theory coupled 
to gravity. We introduce positively and negatively dressed tachyon operators, 

_ r(l + |g|) r , 2 iqX+(--2+\q\)4> y - 1 1! ~ ^ / d 2 i 9 X + (-2-|g|)0 (^OK) 

Vg ~ T(l~\q\)J ' Vg ~ T(l + \q\)J ' ^ 

normalized to remove the usual external leg factors from the correlation functions. For the 
two-, three-, and four-point functions of positively dressed operators in the conventional 
Liouville theory we have [] 

A - —f\li\ r - n .\t-^H\(ll\ + \l2\ + \q3\)/2 

|y»l i 

<t«* = (l - - - • (5.26) 

We may now modify the theory as in example (C), by changing the dressing of the 
operator with momentum qi + q 2 . We assume that this flips the sign of the Liouville energy 
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corresponding to this state appearing in the intermediate channel, and the resulting four- 
point function is 



^™ 4 = ^(E^r 2+(kil+k2l+k3l+k4l)/2 (1 + 



<?i + <72| 9i + 93 ?i + 94 



(5.27) 



Similarly, the change of dressing of the q\ + q% operator changes the sign of the corresponding 
term in the four-point function, etc. It is easy to check, though, that the quantity X^m, 
defined in ( |5.24| ), is invariant under these changes. In fact, 

X qmq 3 q4 = I 9i 1 1 92 1 1 9s 1 1 94 1 • (5.28) 

A change of dressing of one of the external operators is implemented by |gj| — > —\qi\, which 
indeed flips the sign of X. These properties of X qiq2q3qi , found with plausible assumptions 
about Liouville theory, are in complete agreement with our calculations in modified ma- 
trix model. All other normalization independent quantities involving three- and four-point 
functions can be built from (|5.10|) and ( |5.24| ). The correlators of modified operators thus 



agree with the correlators of negatively dressed operators, up to possible rescalings of the 
operators. We expect that this agreement extends to higher-point functions. 

It is interesting that behind the four-point function, which is quite complex, we have 
uncovered a more fundamental object, X qiq2q3q4 , which transforms very simply under changes 
of the Liouville dressing. It would be interesting to check if such simpler objects can be 
defined for higher-point functions. Construction of such objects is reminiscent of the work 
of Di Francesco and Kutasov [] who built the c = 1 correlators in terms of more elementary 
"vertices". Perhaps such objects hold the key to a better understanding of the Liouville 
calculations. 



6 Conclusion 

The recently improved understanding of the modified matrix models with fine-tuned worm- 
hole weights opens the possibility of many new insights into random surfaces. In this paper 
we have calculated some of the simplest modified correlation functions, and there are many 
possible generalizations of our work. It is remarkable that our calculations, which on a sphere 
reduce to Legendre transforms, have a hidden relation to operators with the negative branch 
of Liouville dressing. We hope that more progress will come from a deeper understanding of 
this effect. 



Acknowledgements 

This work was supported in part by DOE grant DE-FG02-91ER40671, NSF grant PHY90- 
21984, the NSF Presidential Young Investigator Award PHY-9157482, James S. McDonnell 
Foundation grant No. 91-48, and an A. P. Sloan Foundation Research Fellowship. C. S. is 
supported by Deutsche Forschungsgemeinschaft. 



19 



Appendix A 



In this appendix we give a careful derivation of the fundamental formula relating the con- 
tinuum one-matrix models and corresponding modified models: 



c 



F(Ut}) = f dt e tt+F(t,{t}) ( A .r 



where {£} = {ti,t 2 , ■ ■ ■} are couplings to scaling operators around a particular critical point. 
Apparently the scaling operators o n conjugate to t n , n > 1 are spectators in the continuum 
formula ( |A.1| ), so that we can extend it to any correlator without punctures 



((Tj-j^ . . . <Ji s )disconn. J dt 6 • • • Cj s )disconn. • (A .2) 



However, it turns out that the matrix model proof of ( |A.1| ) involves some subtleties, due 
to the fact that the discrete scaling operators in the old and the modified models are not 
exactly the same. 

Let us consider a perturbation of the fcth multicritical modified matrix model of the form 
(we use even potentials for simplicity): 

Z k [X,g,{f}] = J Me -<n(.W^^^v.) (A3) 
The multicritical potential is identical to the one of the conventional model, 

^) = E c p $2p > (A.4) 

and A is the bare cosmo logical constant. The bare scaling operators have the general form 

a n = 52gV>Tr&* . (A.5) 

p>i 

The coefficients can be explicitly computed in this model, because the planar limit is 
solvable. Let us review the main points involved in the solution (see ref. []). 

All planar properties can be extracted from the eigenvalue density. In the one-cut phase 
it is given by 

P( x ) = W- E (rn + k + l)A k (R)~g m+k+1 x m VR^ , (A.6) 

m,fc>0 

with the definitions 

g P = c p + f ■ g p - g(-^Tr$ 4 )8 Pt2 , (A. 7) 

MR) = P) (f)' . (A,) 

Since g p depends on the Tr $ 4 condensate we are led to a self-consistent problem given by 
the equation 

— Tr $ 4 ) = J dx p{x) x 4 . (A.9) 
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Fortunately, for the model at hand the condensate enters only linearly in this equation, and 
we can solve explicitly for the feedback. The final answer for the cosmological constant as a 
function of the eigenvalue endpoint is 

^)--i|l + ^-EMp(l-^2)(cp-c 2 +f.?,)| . (A.10) 




1A\ 

Also, the string susceptibility is given by 

y=A/l Tr$ 4\ = A l , An) 

X d\\N / l-gA 2 2 1 J 

For f = a critical point with positive susceptibility exponent 7fc = occurs at R c when 
X~(R-RcT l and \{R) ~ A(i2c) + C(R - R c ) k+l . (A.12) 

If we set R c = 8, then the critical value of g is g c = 1/576. Scaling operators are 
defined as deformations of eq. ( |A.12j ) which do not shift the location of the critical point 



g c and X(R C ). This is an important requirement because, in the end we want the couplings 
f to act as sources, and non-universal dependence on f through the value of the critical 
point invalidates the scaling property. As a result, scaling operators (already in conventional 
matrix models) involve one more tuning than the multicritical potentials. 

From eq. (|A.10|) we find the critical conditions for the coefficients of a n , n > 1: 

X n(R) = ^| Yj> A P (l - J^ A ^) 9p ] ~ ( R ~ 8 ) n+a . (A.13) 

where a is a positive integer. In the conventional models g = and a = 0. Remarkably, 
a = 1 in the modified models. Indeed, the derivative of ( |A.13| ) is 

dXn _i-gA 2 2 v - , _ s . (n) , _ ^„, +Q _i 



dR 2A 2 R 



J2P(P - 2)A p g^ ~(R- S)"^- 1 . (A.14) 
p>i 



At the critical points with positive 7 we have g = 1/576, 1 — gA\ ~ R — 8 and X n (R) is at 
least quadratic in R — 8 for n > 1. This shift is ultimately responsible for the absence of a 
scaling operator with dimension 7fc- From ( |A.14j ) we also see that the scaling operators in 
the modified model differ from those at g — only in the coefficient of Tr $ 4 which does not 
enter ([A. 14 ), and is determined from ([A. 13 ) by requiring stability of the critical point. The 



final result is 



a n = a n + 5gi n) Tr^ , (A. 15) 



where a n is the conventional bare scaling operator, and 



with a spectrum of gravitational dimensions at the kth critical point: 



d(a n ) = (n + l)% = ^±±. (A.17) 
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We see that the bare operators in the two phases are different by a shift of the Tr $ 4 
term. Alternatively, we can work with the same deformations and a shifted cosmological 
constant A — > A + f • 5g 2 . It turns out that the required redefinition of A ensures the delicate 
balance of non-universal terms needed for the scaling of formula (|A.l ) . 

Setting f n = r n we write 

r ■a = T-a + T-5g 2 Ti<5> 4 (A.18) 
and, following ref. [] we introduce the Gaussian representation 

Z k [g,X,{r}} = -JL J + ~dxe-€( x +^-ci-*f Zk[g = , c 2 - x , {r}] . (A.19) 

Note that Z k in the integrand is the partition function of the k-th multicritical model, 
deformed by the conventional scaling operators. Next we separate the non-universal (analytic 
in x) terms on the sphere 

\ogZ k [g = 0, c 2 - x, {r}] = N 2 (-a{ r} x + ^a { 2 } x 2 ) + F(x, {r}, N 2 ) . (A.20) 

The t dependence of the coefficients ai,a 2 can be extracted from the previously studied 
planar solution 

»_ d /! rp..^. OYirt. 



ai -dA(iv Tr *) (9 = ' A = C2) = W If J = 576 - (A ' 21) 

Note that a 2 = 576 independently of r n . On the other hand, a\ depends linearly on r: 

a[ T} = (^Tr $ 4 )(^ = 0, A = c 2 ) = a x (r = 0) + 576 r • 5g 2 = ai + a 2 T-5g 2 . (A.22) 
Finally, we can write the Laplace transform, ready for scaling 

Z k [g, A, {T}] = COnst. X dx e ^a 2 xA+F( X ,{r},N^ ^ ^ ^ 

J — oo 



where A is given by 



M 

\-\ c = \-C2 + T-5g 2 - — = \-C2- — , (A.24) 

a 2 a 2 



and the scaling variable A is independent of the scaling deformations r n . As a result, we 
can safely take derivatives in r n to obtain a discrete version of ( |A.2[ ). We see that the 



slight difference between a n and a n ensures the stability of the critical point A c under such 
deformations. Formula ( |A.1| ) follows then from the scaling 



2 _ 2 2711 

X^tN'i-z, A~tiV^ , Tn^tnN-r-z (A.25) 



With tt- + TT-r 

2—7 2—7 
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